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ABSTRACT 

This paper presents the case of CoRoT LRa02_E2_0121, which was initially classified as a Neptune-size transiting-planet candidate on a relatively 
wide orbit of 36.3 days. Follow-up observations were performed with UVES, Sandiford, SOPHIE and HARPS. These observations revealed a 
faint companion in the spectra. To find the true nature of the system we derived the radial velocities of the faint companion using TODMOR — 
a two-dimensional correlation technique, applied to the SOPHIE spectra. Modeling the lightcurve with EBAS we discovered a secondary eclipse 
with a depth of ~ 0.07%, indicating a diluted eclipsing binary. Combined MCMC modeling of the lightcurve and the radial velocities suggested 
that CoRoT LRa02_E2_0121 is a hierarchical triple system with an evolved G-type primary and an A-type:F-type grazing eclipsing binary. Such 
triple systems are difficult to discover. 

Key words. Planetary systems - binaries: eclipsing - Techniques: photometric - Techniques: radial velocities 



1. Introduction 



Table 1. Coordinates and magnitudes of C0121. 



Since its launch, th e CoRoT space mission (Bagljnetalj|2006; 
lAigrain et alJ l2008t lAuvergne et al.l |2009) has obtained more 
than 10 5 lightcurves of stars in more than 15 different fields that 
are spatially projected toward the center and anticenter of th e 
Galaxy, in search of transiting exoplanets (iDeleuil et al.ll201 Th . 
In order to confirm a new transi ting planet, a candidate has to 
pass a sequence of tests (e.g., Brown 200 3b. The tests i nclude 
careful analysis of the lightcurv e dCarpano et al.1 12009). pho- 
tometric follow-up observations ( IDeeg et alJl2009h, and f inally 
spectroscopic follow-up observations (Bouchy et al. 2009). 

The main source of false positive alarms in the CoRoT 
sample are eclip sing binary systems in various configurations 
(lAlmenara et al.l 1200 91. One possible strong indication of the 
binary nature of a system is the detection of a faint compo- 
nent in its observed spectrum. However, sometimes the stellar 
secondary is too faint to be noticed at first glance, and spe- 
cial t ools are needed to detect it. Such a tool is TO DMOR 
(e.g., IZucker&Mazehlll994t IZucker et al.ll2003l 12004b . a two- 
dimensional correlation algorithm that was proved several times 
in the past to be an efficient tool for measuring the radial veloc- 
ities (=RVs) of the two components of a binary system, even 
with a relatively f aint secondary (e.g., Maz eh et aDl 1995lll997t 
iTorres et al.lfl995h . 

This paper presents the case of the transiting Neptune-size 
candidate LRa02_E2_0121 (= C0121). By careful analysis of the 
CoRoT lightcurve combined with spectroscopic follow-up ob- 



Main identifiers 


CoRoT ID 


110680825 




2MASS ID 


06520760-0526137 




Coordinates 


RA (J2000) 


06"52 m 07 J .61 




Dec (12000) 


-05°26'13".8 




Filter 


Magnitude 


Source 


B 


12.32 ± 0.05 


ExoDat" 


V 


11.75 ±0.08 


ExoDat 


r' 


11.37 ±0.06 


ExoDat 


V 


10.95 ± 0.07 


ExoDat 


.1 


9.723 ± 0.025 


2MASS 6 


H 


9.197 ±0.020 


2MASS 


K 


9.102 ±0.021 


2MASS 



ueieuii ei ai. izui 
' Cutri et al. (2003 



* Based on observations made with the 1.93-m telescope at 
Observatoire de Haute-Provence (CNRS), France, the 3.6-m telescope 
at La Silla Observatory (ESO), Chile (program 184.C-0639), the VLT 
at Paranal Observatory (ESO), Chile (program 083.C-0690), and the 
2. 1-m Otto Struve telescope at McDonald Observatory, Texas, USA. 



servations of this star, we concluded that it is a hierarchical triple 
system with an evolved G-type primary and a grazing eclipsing 
binary. The case is presented here to illustrate the long curvy 
road one might need to travel until the true nature of a system is 
revealed and to demonstrate the potential of TODMOR, which 
can reveal the diluted-binary nature of a system with a relatively 
modest investment of observational resources. 

Section 2 gives some details about the star and presents the 
CoRoT white lightcurve. Section 3 describes the spectroscopic 
follow-up observations and the spectral analysis. Section 4 de- 
scribes our application of TODMOR to the observed SOPHIE 
spectra. Section 5 presents our analysis of the CoRoT lightcurve 
and the combined MCMC modeling of the photometric and RV 
data. The astrophysics of C0121 as a hierarchical triple system 



1 



Tal-Or et al,: CoRoT LRa02_E2_0121: Hierarchical triple system with a giant primary 



-15 



o) .10 



£ 

> 
jo 



10 - 



1 1 1 1 1 1 1 1 1 1 1 



263.4 263.6 263 8 



299.6 299.8 300 300 2 



3358 336 336 2 336.4 



250 260 270 280 290 300 310 320 330 340 350 
HJD- 2454545 

Fig. 1. White lightcurve of C0121. Top: the calibrated lightcurve 
after rebinning oversampled points back to 512 s and removing 
outliers. The calculated long-term trend is overplotted with a 
solid black line. Middle: the detrended lightcurve (for conve- 
nience the calibrated lightcurve was moved up by lOmmag). 
Insets: zoom on the three individual transit-like events. 



with a G-type giant is detailed in Section 6. Section 7 presents 
some general conclusions derived from this specific case. 

2. The CoRoT lightcurve 

The star C0121 was disc overed as a transiti ng candidate by the 
alarm mode of CoRoT ( Sura ce et al. 2008) during the LRa02 
run, which lasted from Nov. 16, 2008 to Mar. 11, 2009. Its 
white lightcurve displayed three transit-like events of about 
0.3% depth with a period of about 36.3 days, and therefore 
C0121 was identified as an interesting candidate of a Neptune- 
size planet with a relatively long period. 

Figure Q] shows the calibrated and detrended white 
lightcurves of this candidate, with a zoom on the three transit- 
like events in the insets. Trends were removed with local regres- 
sion that fits a second-degree polynomial to one-day-long data 
with weighted linear least squares. Table Q] lists some basic in- 
formation on C0121. 



3. Spectroscopic follow-up observations 

3.1. Radial velocity measurements 

The first five RV measurem ents of C0121 were obtained with 
the SOPHIE spectrog raph dBouchv & The Sophie Teamll2006t 
iPerruchot et alj l2~008h . mounted at the 1.93 m telescope at the 
Observatoire de Haute Provence, France, between Feb. 26 and 
Apr. 04, 2009. Following the discovery of the transits in the 
alarm mode, the observations began before the end of the LRa02 
run and spanned 37 days to cover a full orbital period of the tran- 
siting object. The observation timings were chosen to be close 
to the quadrature times of the transiting object to maximize the 
chance that small RV variability will be detected. These five ob- 
servations gave us the first signs of the complexity of the system. 
Instead of the expected 36.3-day periodicity, we detected a long- 
term drift, with an amplitude of about 0.5 km s _1 . 

Six more observations of C0121 were obtained when the 
field was observable again, between Nov. 05 and Dec. 04, 2009. 
Three were acquired with SOPHIE and three other ones with 




200 250 
HJD - 2454800 



Fig. 2. RVs of the primary in C0121 as measured by the SOPHIE 
(black circles) and HARPS (red boxes) spectrographs using 
CCF The solid black line is a second-degree polynomial fit. The 
residuals are plotted at the bottom. 



the HARPS spectrograph dPepe et alJl2000t iMavor et al.ll2003l) . 
mounted at the ESO 3.6 m telescope in La Silla Observatory, 
Chile. All spectra had a typical signal-to-noise ratio (S/N) of SO- 
SO per pixel at 550 nm. In the additional six spectra we found 
the same long-term drift, which at that time had an observed 
amplitude of about 7 kms _I . The long-term variation indicated 
that C0121 is a member of a binary system with a long pe- 
riod, but did not rule out the existence of a planet in the system. 
There could still be a planet orbiting one of the members of the 

J II 

wide binary, as in the cases of 1 6 Cyg B (Cochran et al. 1997), 
y Cep hei A (lHatzes et al.l 120031) . or 30 Ari B dGuenther et al.1 
2009). Obviously, one needs more observations to understand 
such a complex system. Figure [2] shows the RVs of the primary, 
indicating that the orbital period is too long for us to derive an 
orbit, so we fit the measured RVs with a second-degree polyno- 
mial. 

The plot started to thicken when we noticed that the three 
spectra taken with HARPS revealed a faint secondary. Following 
this, we re-analyzed the eight spectra already obtained with 
SOPHIE, fitting a sum of two Gaussia ns to the cross correlation 
function (= CCF, Baran neet al.llT9 96) derived with a G2 mask 
dPepe et al.ll2002l) . Due to the low S/N of the secondary feature 
in the CCF, we succeeded in measuring its RV in only 9 out 
of the 1 1 measured spectra — three from HARPS and six from 
SOPHIE spectra. All measured RVs are listed in Table [2] 

The nine measured velocities of the faint secondary were 
consistent with the 36.3-day period found by CoRoT, but with 
a large RV semiamplitude of ~ 62 km s _1 . This finding ruled out 
completely the existence of a transiting planet in C0121, and in- 
dicated that the extremely shallow transit-like signal is coming 
from its faint companion, which is by itself a diluted eclipsing 
binary. 

Along the way we found one more surprise: the measured 
RVs of the faint companion of C0121 were in anti-phase with 
the CoRoT lightcurve; i.e., the detected eclipses occurred on the 
rising part of the RV curve. This means that the eclipses detected 
in the CoRoT lightcurve are caused by the star whose spectrum 
is detected in the composite spectra of the system, passing in 
front of its companion. The spectrum of the eclipsed star is not 
seen in the observed spectra. This is a very unusual case, because 
in most eclipsing binaries the primary eclipse occurs when the 
primary star, whose spectrum dominates the spectra, is eclipsed. 

It was clear then that C0121 includes at least three stars. To 
facilitate the discussion we denoted the primary star in the spec- 
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Table 2. The measured RVs of C0121 (kms 1 ). 



HJD 


CCF: 


CCF: 


S/N 


Per Pixel 


(-2450000) 


Primary 


Secondary 


(at 


550 nm) 


4889.30773 


30.091 +0.034 


70.07 ± 


1.78 




33.9 


4904.29007 


29.789 ± 0.035 








35.1 


4911.38470 


29.694 ± 0.045 


-45.34 ± 1.76 




27.8 


4914.30337 


29.663 ± 0.023 


-21.05 i 


= 1.34 




49.8 


4926.31911 


29.403 ± 0.030 


74.39 ± 


2.05 




43.3 


5140.63895 


24.128 ±0.042 








29.4 


5141.70360 


24.071 ±0.039 


72.67 ± 


2.52 




31.3 


5148.56025 


23.815 ±0.029 


62.39 ± 


2.15 




37.8 


5164.81355 


23.312 ±0.012 


-41.98 ±0.87 




51.0 


5168.84324 


23.225 ± 0.014 


-11.56 ±0.74 




48.4 


5169.70123 


23.166 ±0.016 


-0.86 ± 


0.73 




43.8 


HJD 


TODMOR: 


TODMOR: 


S/N 


Per Pixel 


(-2450000) 


Primary 


Secondary 


(at 


550 nm) 


4soy. JU / 15 


3U. 11 ± U.l 1 


74.08 ± 


1.32 




to n 
ii.y 


4904.29007 


29.82 ±0.12 


-27.05 4 


= 1.14 




35.1 


4911.38470 


29.72 ±0.19 


-44.60 4 


= 1.77 




27.8 


4914.30337 


29.68 ±0.10 


-22.13 ±0.89 




49.8 


4926.31911 


29.40 ±0.12 


75.61 ± 


1.20 




43.3 


5140.63895 


23.93 ±0.11 


68.04 ± 


1.33 




29.4 


5141.70360 


23.85 ±0.11 


73.55 ± 


1.28 




31.3 


5148.56025 


23.89 ±0.10 


65.67 ± 


1.20 




37.8 



tra by C0121-A, the secondary star in the spectra, whose RV 
was observed to vary with a large amplitude by C0121-Bb, and 
the eclipsed star during the primary eclipses detected so far by 
C0121-Ba. 



3.2. High-resolution spectroscopic observations 

A reconnaissance high-resolution (R « 47 000) spectrum 
of C0121 was acquired on February 13, 2 009, using the 
Sandiford Cassegrain echelle spectrograph (iMcCarthv et al.l 
1 1993b mounted at the 2.1 m Otto Struve Telescope at McDonald 
Observatory. The analysis of the spectrum revealed a G-type 
giant star with r eff = 5500 + 250 K, logg = 3.0 ± 0.2, and 
vsin/ =12 + 1 kms -1 . 

Furthermore, high-resolution, high-S/N spectroscopy 
of C0121 was performed with the UVES spectrograph 
dDekker & D'Odoricd [19921) mounted at the ESO 8.2 m Very 
Large Telescope (= VLT) in Paranal Observatory, Chile. Two 
consecutive exposures of 3000 sec each were acquired on April 
8, 2009, under the ESO program 083.C-0690. The adopted 
set-up yielded a resolving power of R as 67 000 in the spectral 
range of 3280 - 6820 A, with a final S/N of about 320 per 
pixel at 5500 A. Following the proced ure already adopted 
in some CoRoT di scovery papers (e.g., iDeleuil et al.l [2008; 
iGandolfiet all 120101) . the UVES spectrum was used to derive 
the fundamental photospheric parameters of C0121. Ignoring 
the additional light from C0121-Ba and C0121-Bb, we found 
r eff = 5640 + 100 K, logg = 2.9 ± 0.1, [M/H] = -1.0 ± 0.2, 
and vsin/ = 11.5 + 1.0 kms -1 , corresponding to a GO III star. 
Allowing a dilution factor and fitting only the Mg I, Fe I/II and 
Cal lines we found 7/ eff = 4980 + 100 K, logg = 2.6 + 0.1, 
[M/H] = -0.5 + 0.2, and vsin/ = 11.0 + 1.0 km s" 1 , 
corresponding to a G8 III star. 

The two observations indicate that C0121-A is probably a 
G-type giant star. 



triple system with a giant primary 




Radial Velocity [km/s] 

Fig. 3. TODMOR plot for the first exposure of CO 1 2 1 . The upper 
and lower panels show the primary and secondary cuts through 
the two-dimensional correlation function. 



4. TODMOR analysis of the SOPHIE spectra 

To better understand the unusual C0121 system, we analyzed 
all eight SOPHIE spectra with our two-dimensional correlation 
algorithm TODMOR, and derived the RVs of C0121-Bb from all 
eight spectra. Our best model included as the primar y template 
the o bserved spectrum of HD 32923, a G4V star (Whit e et all 
120071) . rotationally broadened by 9 kms 1 , and as the secondary 
templ ate the observed s pectrum of HD 185395 (9 Cyg), an F4V 
star dKupka et al.l 12004). The secondary/primary flux ratio was 
~ 0.07. ~ 

Despite our efforts, we did not detect any spectral signa- 
ture of C0121-Ba, even when analyzing the data with TRIMOR 
(Maz eh et al.ll2009l) . an algorithm to analyze multi-order spectra 
of triple systems. This might indicate that Ba is a substantially 
hotter star, with spectral type A or even earlier. Such stars have 
only few broad spectral lines in the optical band, and therefore 
it is hard to measure their RVs, especially when their light is 
blended with the light of other stars. 

Figure [3] shows the TODMOR results for the first SOPHIE 
exposure of C0121. The upper and lower panels of the figure 
show primary and secondary cuts through the two-dimensional 
correlation function that run through the correlation peak. The 
primary cut (upper panel) is parallel to the primary RV axis, 
freezing the secondary velocity at its derived velocity. The sec- 
ondary cut (lower panel) runs parallel to the secondary RV axis, 
freezing the primary velocity at its derived velocity. 

We note that in the upper panel the correlation drops by ~ 0.6 
when moving away from the peak, because we change the veloc- 
ity of the primary template in the model. On the other hand, the 
correlation in the lower panel drops only by ~ 0.002, because we 
change the velocity of the secondary template, which contributes 
only ~ 7% of the light, and the primary velocity is kept at its best 
value. Nevertheless, the peak at the lower panel is prominent, in- 
dicating the high significance of the detection of the secondary 
velocity. 

The TODMOR RVs of C0121-A and C0121-Bb are listed in 
Table [2] Figure presents the RVs of C0121-Bb with our best 
orbital solution, which is described in the next section. 

We also estimated the effective temperatures of C0121-A and 
C0121-Bb from the SOPHIE spectra by try ing a range of high- 
resolu tion PHOENIX synthetic models ( Hau schildt & Baronl 
1 19991) with different T eB values as templates. TODMOR found 
T'eff.A = 5650 + 200 K and r eff ,B b = 6500 + 500 K, in agreement 
with the analysis of the UVES and Sandiford observations. 
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Fig. 4. RVs of C0121-Bb as measured by the SOPHIE (black 
circles) and HARPS (red boxes) spectrographs. SOPHIE RVs 
were derived with TODMOR, while HARPS RVs were derived 
by fitting a sum of two Gaussians to the CCF. The solid line 
is the Keplerian model produced from the orbital parameters in 
Table [3] and the dashed line is the center-of-mass velocity. The 
residuals are plotted at the bottom. 



5. Modeling the eclipsing binary C0121-B 

To estimate the photometric parameters of C0121-B we 
re-analyzed the CoRoT white lightcurve with EBAS 
dTamuz et all l2006t iMazeh et all 120061) , an algorithm for 
analyzing eclipsing-binary ligh tcurve based on the EBOP sub- 
routines (Popper & Etzell ll98ll) . Our EBAS analysis discovered 
a faint secondary eclipse with a depth of (0.07 + 0.01)% at the 
orbital phase ~ 0.52, showing again that C0121-B is actually an 
eclipsing binary (see Figure 01. The CoRoT lightcurve covered 
three primary and four secondary eclipses. 

Final parameters were derive d with MCMC analysis (e.g., 
iTegmarket ail l2004t iFordl 120051) of the white lightcurve, to- 
gether with the RV data. The input data were the detrended white 
lightcurve from CoRoT, which included 19,348 points and the 
11 measured RVs of C0121-Bb. From the RVs listed in Table 
|2] we used the 8 TODMOR-derived SOPHIE RVs and the 3 
CCF-derived HARPS RVs. Errors for the photometric data were 
derived from the local scatter of the detrended lightcurve (see 
Figure [TJ. 

The sizes of the MCMC Gaussian perturbations were set 
by the error estimations of EBAS and by a relatively short 
MCMC run (of ~ 10 5 accepted steps). Following the Metropolis- 
Hastings algorithm, trial points with lower x 1 were accepted, 
whereas trial points with higher x 1 were accepted only with 
probability of e xp(-Ay 2 /2 ), assuming the observational errors 
to be Gaussian (Ford 2005). The^ 2 for each trial point was the 
sum of the lightcurve^ 2 and the^ 2 of the RV data. We used flat 
prior distributions for all free parameters of the model, all de- 
tailed in Table 3. The absolute values of e cos a>, e sin to, J s , r t , k, 
x and Lj were constrained to be between and 1 . 

The fractional sum of radii (r,) and the impact parameter (x) 
can be expressed by the more conventional parameters as 



R„ + R s 



JS 10' 

o 



-0.01 -0.005 



0.005 0.01 0.51 0.515 0.52 0.525 0.53 

Phase 



Fig. 5. Zoom on the primary and secondary eclipses in the phase- 
folded lightcurve of C0121 with our best model overplotted with 
a solid red line. The residuals were moved down by 10 mmag for 
convenience. 



and 



1 



r, 1 + e sin d 



(2) 



(1) 



where R p and R s are the primary and secondary radii, a is the 
orbital semimajor axis, i the i nclination, e the e ccentricity, and a> 
the longitude of periastron dTamuz et al. 2006). Using this con- 
vention, x = 1 when the components are grazing but not yet 
eclipsing. We used linear lim b darkening law, with coefficients, 
u p and u s , fixed according to Singj d201C ), assuming an A5 and 
F5 dwarfs for the primary and secondary stars in the eclipsing 
binary. 

The MCMC run consisted of 2 10 6 accepted steps. The de- 
rived value of each parameter was its MCMC median, and its 
confidence limits were estimated as the range of MCMC values 
that cover the central 68.3% of the chain. 

Table [3] gives the full list of the estimated parameters of 
C0121-B from the MCMC analysis. Most of the model pa- 
rameters showed normal posterior probability distribution. For 
the ones with skewed distribution, Table [3] gives uneven con- 
fidence limits. Figure [5] presents the folded lightcurve and our 
best model, zoomed on the primary and secondary eclipses. 

We note that the high x value we found indicates that the 
C0121-B eclipsing binary is grazing. In addition, the signifi- 
cantly non-zero center-of-mass acceleration value we found in- 
dicates that the C0121-B binary is itself a member of a hierar- 
chical triple system. 

The time of periastron passage (T), as well as some other pa- 
rameters of C0121-B, which are listed at the bottom of Table[3] 
were analytically calculated for each point in the MCMC chain. 
Their values and errors were estimated the same way as for the 
free parameters of the model. Only for the Bb:Ba flux ratio did 
we take the most probable value of J s k 2 as our estimation, where 
J s is the Bb:Ba surface-brightness ratio and k their ratio of radii 
(see TableO, and not the median one, because we found its prob- 
ability distribution to be extremely skewed. 

C0121 was bright enough for CoRoT to observe it with 
its th ree different bands (e.g., lAigrain et alJl2008t iBorde etaD 
l2010h . The reduction of the colored lightcurves, which was done 
after the completion of the LRa02 run, revealed different depths 
of the primary eclipse — (0.42 ± 0.02)% in blue, (0.35 + 0.02)% 
in green, and (0.18 + 0.02)% in red, indicating again that C0121- 
Ba is probably hotter than the other two stars in the system. 
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Table 3. The parameters of C0121-B from the combined lightcurve and RV modeling. 



Symbol 


Parameter name 


Value 


Units 


Orbital parameters 


P 


Period 


36.2873 ± 0.0025 


day 


ecosoj 


Eccentricity x Cosine Longitude of periastron 


0.0331 ±0.0003 




es'ma) 

K 

7 


Eccentricity x Sine Longitude of periastron 

RV semiamplitude 

Center-of-mass RV at JD = 2454889 


029+o.oos 

Q007 

62.2 ±0.5 
14.5 ±0.6 


km s" 1 
kms -1 


7 


Center-of-mass acceleration 


20.6 ± 2.8 


ms _1 d _I 


T 


Time of periastron (calculated analytically) 


2454876.7 ± 0.7 


JD 


Photometric parameters 


To 


Time of center of primary eclipse 


2454844.9152 ±0.0020 


JD 


J* 

n 


Surface-brightness ratio (secondary/primary) 
Fractional sum of radii 


0.61 ± 0.23 
0.049 ± 0.006 




k 

X 


Ratio of radii (secondary/primary) 
Impact parameter 
Third-light (blending) fraction 


0.60 ± 0.25 
0.88^ 
0.86 ± 0.04 




Up 


Limb-darkening coefficient of primary (fixed) 


0.55 






Limb-darkening coefficient of secondary (fixed) 


0.58 




Other estimated parameters 


i 


inclination (see Eq.(2)) 


87.45-- 


degree 


f 
1 


Semimajor axis of secondary 
Mass function of Bb 

Mass ratio (M Ba /M Bb ) assuming M Bb = 1.3M S 


0.2076 ±0.0017 
0.903 ± 0.022 
1.73 ±0.03 


AU 

M G 


FBb'-FBa 


The Bb:Ba flux ratio 


6+ 03 





6. The C0121 system 

The light of C0121 is coming from (at least) three stellar sources: 
a G-type star (C0121-A) on a wide orbit with an unknown period 
and an eclipsing binary (C0121-B) that includes an F-type star 
(C0121-Bb) orbiting a hotter and more massive star (C0121-Ba) 
on a nearly circular orbit of ~ 36.3 days. In addition, the center 
of mass of the C0121-B system is itself on a wide orbit with un- 
known period. We explore now the assumption that the three ob- 
served stars, A, Ba, and Bb, consist of a single hierarchical triple 
system. To do that, we use the following previously-estimated 
quantities and relations, derived from the white lightcurve and 
spectra of C0121: 

- The effective temperature of C0121-Bb is 6500 + 500; 

- The effective temperature of C0121-A is 5650 + 200; 

- The Bb:Ba surface-brightness ratio is 0.61 + 0.23; 

- The Bb:Ba flux ratio is 0.6+°- 3 ; 

- C0121-A contributes 86 + 4% of the total light of the system; 

- The Ba:Bb mass ratio q depends on Mbi, through 
{Msbf 1 sin 3 i) q 3 - q 2 - 2q - 1 =0, where the mass function 
/ is 0.903 + 0.022 M and the inclination i is 87.4 + 0.4°; 

- The center-of-mass acceleration of C0121-B is 20.6 + 
2.8ms- 1 d- 1 ; 

- The average acceleration value of C0121-A is -25 + 

lms-'r 1 . 

To proceed, we assume that Bb is on the main sequence. From 
Y 2 st ellar isochrones of solar-metallicity stars (Demaraue et al] 
2004) we take for it the following values (see Fig. [6]): 

- The My brightness of Bb is 3.7; 

- The mass of Bb is 1.3 M . 

To estimate r e ff Ba we used the derived surface brightness ra- 
tio and calculated the expected surface-brightness ratio between 
various PHOENIX synthetic models, taking the CoRoT total 
efficiency (Mi chel et alj l2QQ9h into account. We got r e ff,Ba = 
7400^ggy° K. From the assumed My value of Bb, and using the 
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12,000 10,000 8,000 6,000 4,000 
T„["K] 



Fig. 6. Y 2 stellar isochrones from Demar que et al] d2004l) for 
solar metallicity ([Fe/H] = 0.05, [a/Fe] = 0) and ages of 
0.2,0.4,0.6,0.8, and 1 Gyr, from top-left to bottom-right. Left: 
Mv/(T e ff) isochrones. Right: M(T e ff ) isochrones. Our estimations 
forC0121-A, C0121-Ba, and C0121-Bb are marked by a red di- 
amond, a green square, and a blue circle, respectively. The abso- 
lute magnitude and mass of C0121-Bb were taken as the average 
value of the isochrones for a 6500 K star, assuming no error. 



derived brightness ratios, we got for C0121-Ba and C0121-A an 
My values of 3.1^^ and Q.6_A, respectively. From the assumed 
mass of Bb and the derived ratios, including the acceleration ra- 
tio between A and B, we got for C0121-Ba and C0121-A masses 
of 2.25 + 0.04 M and 2.8 + 0.4 M , respectively. 

To see whether the three temperatures, brightnesses, and 
masses can be fitted on a single isochrone, we plotted the Y 2 
stellar isochrones for solar metallicity and ages of 0.2-1 Gyr in 
Figure [6] together with our estimations for the effective temper- 
atures, masses, and My values of the three observed stars. It can 
be seen that a single isochrone of ~ 0.6 Gyr could fit all nine 
estimated values, supporting the hierarchical triple hypothesis. 

Another way to estimate the stellar temperatures is to use the 
lightcurves of the different CoRoT bands. We derived the six J s j 
and Lij values, where j indicates the color, using MCMC as be- 
fore, and searched for the best temperatures of the three stars, 
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together with the best Ba:A ratio of radii, fc^A, similar to the 
way we estimated r e ff,B a from the white lightcurve parameters. 
To do that we had to estimate the wavelength range of the blue, 
green, and red CoRoT bandpasses. From the observed magni- 
tudes of C0121 in each channel, we estimated the two effective 
wavelengths that separate the channels in the CoRoT mask to 
be roughly at 5100 and 5600 A. We got T eSA = 5700 + 200 K, 

Trft.Ba = 7200 ± 400 K ' r eff3b = 6500 ± 300 K, and k BaA = 
0.20 + 0.02. These values agree with the above estimated values 
and confirm that C0121-A is indeed a G-type giant star. 

The visual and absolute magnitudes of C0121-A yield a dis- 
tance of 1 .7*H kpc to the system, consis tent with being in the 
Perseus arm of the Galaxy (Russeil 2003b. 

7. Conclusions 

The available photometric and spectral information for C0121 
indicates that the system is a hierarchical triple system and 
that the distant companion, C0121-A, is a G-type giant star. 
Hierarchical triple systems with the distant component in its 
giant phase are quite difficult to discover. Even if the close 
pair is an eclipsing binary, the dilution factor makes the eclipse 
rather shallow. In our case the (graz i ng) ec lipse depth is only 
~ 0.3%. Very recently. iDerekas et alJ (1201 li) discovered a simi- 
lar case, HD 181068, for which the observed depth of the close- 
pair eclipse was also much less than 1%. Not surprisingly, the 
triple nature of HD 181068 was also discovered by a lightcurve 
obtained from space, this time with the Kepler mission. Ground- 
based observations could not have detected the eclipses of the 
close binary in these systems. When the close pair is not eclips- 
ing it is almost impossible to detect the triple nature of the sys- 
tem if the distant companion is a giant, even if the wide orbital 
period is short enough to be discovered by radial-velocity obser- 
vations of the giant. The discoveries of these two systems indi- 
cate that triple sy stems might not b e rare, as al ready found by 
iMavor & Mazehl (ll987h . [Tokovininl (12008b . and lRaghavan etalJ 
(I2010h . 

The present case demonstrates how a promising Neptune- 
size transiting-planet candidate turned out to be a grazing eclips- 
ing binary in a hierarchical triple system. Indeed, hierarchical 
triple systems are expected to be the main source for false pos- 
itive detections of N eptune-size transiting-pla net candidates in 
the Kepler data too (iMorton & Johnsonll201 ll) . The true nature 
of the system could have been revealed by analyzing the first few 
SOPHIE spectra with TODMOR, suggesting that an automated 
TODMOR pipeline could be useful in the long planet-hunting 
procedure. We hope to construct such a pipeline in the near fu- 
ture, so we can detect disguised hierarchical triple systems as 
early as possible. 
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